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Two-dimensional Bloch electrons in a uniform magnetic field exhibit complex energy spectrum. 
When static electric and magnetic modulations with a checkerboard pattern are superimposed on 
the uniform magnetic field, more structures and symmetries of the spectra are found, due to the 
additional adjustable parameters from the modulations. We give a comprehensive report on these 
new symmetries. We have also found an electric-modulation induced energy gap, whose magnitude 
is independent of the strength of either the uniform or the modulated magnetic field. This study 
is applicable to experimentally accessible systems and is related to the investigations on frustrated 
ant if erromagnet ism . 

PACS numbers: 72. 20. My; 73.21.-b; quantum wells, mesoscopic, and nanoscale systems 73.43.-f; 75.10.Jm 



I. INTRODUCTION 

When the spectrum of a two-dimensional (2D) Bloch 
electron in a uniform magnetic field is plotted in the 
energy-flux diagram, a self-similar structure with fractal 
property emerges^ Such a complex structure, called the 
Hofstadter spectrum, arises due to the commensurabil- 
ity between two length scales in this system: the lattice 
constant and the cyclotron radius. The Hofstadter spec- 
trum is one of the earliest predictions of fractal struc- 
ture in solid-state physics. Subsequently, it was found 
that not only the energy spectrum has self-similarity, the 
wave function also exhibits scaling behavior and can be 
analyzed using the renormalization group. 2 

Because of its beautiful structure, the Hofstadter spec- 
trum has attracted many researchers' attention, and the 
spectra for different 2D lattice symmetries have been re- 
ported. Besides the square lattice, there are also triangu- 
lar lattice^ honeycomb lattice, 4 Kagome lattice, 5 and a 
bipartite periodic structure with hexagonal symmetry,^ 
These are all studied within the framework of Bloch 
electrons in a uniform magnetic field, usually assum- 
ing nearest-neighbor (NN) couplings t% only. Including 
and varying the next-nearest-neighbor (NNN) couplings 
ti leads to band-crossings accompanied by exchange of 
quantized Hall conductances between bands For a 
square lattice, detailed scaling analysis reveals a bicritical 
point at t\ = 2*2) accompanied by interesting topologi- 
cal change of the spectrum^ Spectra for systems with 
couplings beyond next-nearest neighbors have also been 
studiedjAS In addition, the external magnetic field, rather 
than being uniform, can be periodically modulated with 
a pattern unrelated to the original lattice. The simplest 
situation when a magnetic lattice overlaps with the elec- 
tric lattice is realized when a ferromagnetic grid is de- 
posited on a semiconductor >ii The interfacial stress be- 
tween two materials would naturally induce an electric 
grid with the same period and symmetry as the ferromag- 
netic grid. More generally, there can also be a magnetic 



modulation with the pattern of a ID strip^i 2 or a 2D 
checkerboarcU^i^ii^ superimposes on the electric square 
lattice. The checkerboard configuration has been real- 
ized experimentally using a superconducting Nb-network 
with periodic magnetic Dy-islandssAS, The calculations of 
the Hofstadter spectra provide the basis to study such 
articifial networks. 

A direct observation of the Hofstadter spectrum has 
been realized using microwaves 1 ^ or acoustic wavesii 
transmitting through an array of macroscopic scatters. 
However, a fractal electronic spectrum is significantly 
more difficult to be realized in an usual solid, whose lat- 
tice constant is only a few angstroms, and a magnetic 
field of the order of 10 4 Teslas is required. For 10 Tes- 
las or less, we can only probe the part of the Hofstadter 
spectrum that reduces to the familiar Landau levels with 
roughly equal spacings in energy. In the last decade, 
different superlattice structures with much larger lattice 
constants are used to cope with this high- field problem. 18 
Besides, several physical systems are closely related to 
the Hofstadter problems and offer alternative angles of in- 
vestigation, for example, the studies of a superconductor 
in a vortex stated a superconducting network in a mag- 
netic field^ and a junction of three quantum wires^S 
Furthermore, recent advance on optical lattices makes 
it possible to implement a lattice Hamiltonian resem- 
bling the effects of magnetic fields with neutral atoms. 21 
This offers great opportunities since not only the mag- 
netic field, but also the lattice symmetry, the potential 
strength, and the relative importance of many-body ef- 
fects can be adjusted in such a system. 

This paper is motivated by a study very different from 
those mentioned above. In a recent paper^ Misguich 
and coworkers, by using the hard-core bosons to represent 
the spin degrees of freedom, and using the Chern-Simons 
transformation to transmute bosons to fermions, mapped 
a 2D frustrated antiferromagnetic problem to a Hofs- 
tadter problem. This approach is subsequently used to 
study the magnetization properties of the J\ — J2 Heisen- 
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berg model on a square latticed After suitable mathe- 
matical mappings and a mean-field approximation, the 
magnetization problem can be reduced to a Hofstadter 
problem with both electric and magnetic super-structures 
superimpose on the original lattice. This motivates us to 
consider the checkerboard super-structure, which is re- 
lated to the Neel phase in the magnetization problem, 
with congruous electric and magnetic modulations. Cou- 
plings up to next nearest-neighbors are considered, which 
are essential to cause magnetic frustration in the J\ — J2 
Heisenberg modeli24 

In this paper, we make a comprehensive survey of the 
symmetries of the Hofstadter spectra with field modu- 
lations. Some of the symmetries already exist without 
modulations, such as the ones related to reversing the 
direction of the magnetic flux, and shifting the flux in a 
plaquette by two flux quanta (see items II and III in Sec. 
II). Some of the other symmetries that are closely re- 
lated to field modulations are reported for the first time. 
In particular, when the system is subject to staggered 
7r-fluxes, its spectrum in the E — cf> diagram has an up- 
dowm symmetry even with NNN couplings (see Fig. 5), 
which is quite unexpected since NNN couplings usually 
destroy such a symmetry^ Besides the studies on sym- 
metries, for systems with only NN couplings, we find a 
simple algebraic connection between the spectra with and 
without the electric checkerboard field. We also find a 
flux-independent energy gap induced by the electric mod- 
ulation (see Fig. 3), which can be explained by using the 
algebraic relation just mentioned. 

This paper is organized as follows. Theoretical formu- 
lation on the system with checkerboard super-structure, 
as well as the discussion of the symmetries of the energy 
spectra, can be found in Sec. II. Major features of the 
Hofstadter spectra are discussed in Sec. III. We summa- 
rize and conclude our results in Sec. IV. The proofs on 
checkerboard-translation symmetry of the spectrum and 
on the existence of the flux-independent energy gap are 
given in the appendices. 

II. THEORETICAL ANALYSIS 

The tight-binding Hamiltonian describing the motion 
of an electron in a magnetic field is given by 

H = -\ E (*« e ** ftu + h - c ) + E Kfih < C 1 ) 

(id) 1 

where 9%j is the magnetic phase factor. For clarification, 
we will replace the label i by (n, m) in the following, 
which denotes the (n, m) plaquette as well as the lattice 
point at the lower left corner of the plaquette. Without 
loss of generality, we take the uniform part of the mag- 
netic flux through plaquette as cj> — 2irp/q with relative 



prime integers p and q. For the checkerboard modulation, 
we have ^/2tt = -A (-l)"+ m and V, = A y (-l)"+ m 
(see Fig. 1). The Landau gauge is used such that the 
magnetic phase factors become 

{^n+l, m; n. m 3 
9 n , m+l; n, m = H + (-l)"+ m 7rA , (2) 
#n+l, m+l; n, m = On, m+l; n+1, 771 = \Tl + 5) 4> ■ 

Due to the modulation in the y direction and under the 
gauge choice in Eq. @ , the tight-binding Hamiltonian in 
Eq. ||TJ now becomes invariant under the y-translation 
m — * m + 2. Thus the Bloch theorem gives 

fn, m — C y C n ^ m(^y) (3) 

for \k y \ < 7r/2, where c n , m+2 {k y ) = c n , m (k y ) and 
Cn,m(ky + 7r) = c n ^ m (k y ). Therefore, the generalized 
Harper equation becomes 

(k y ) + B n c n+1 (k y ) + B n _ic^_i(fcj,) 
= Ec n (k y ) (4) 

where c n (k y ) = (c„, i(k y ),Cn, 2(k y )) T and 

/ -(-l)»Ay ~h cos( Xn )e iS " \ , , 
n ~ {-t lC o S (xn)e- iS " (-l)"Ay J W 

tj ( -t2COs(77„) \ . , 

n ~ \-t2COs( Vn ) -tl/2 ) W 

with \n = ncj) + k y , S n = (— l)™7rA<£, and rj n = (n + 
1/2)0 + k y . It can be easily checked that 

A„ + q = An, B„ + q = B„, (7) 

where Q = q (2q) for an even (odd) integer q. Thus, the 
n in Eq. Q satisfies the condition 1 < n < Q. Besides, 
because of the magnetic translation symmetry, the prim- 
itive unit cell is consisted of q (Q) plaquettes without 
(with) checkerboard modulation. 

The Bloch condition along the x direction can be writ- 
ten as 

Cn, m(ky^) — C V'ti, m(k x , ky} (8) 

for \k x \ < ir/Q, where ^ n +Q, m(k x , k y ) = ip 

71 , 777 ( k X ? ky ) 

and ip n> m (k x + 2n/Q, k y + n) = ip n , m (k x ,k y ). Now we 
only need to solve the problem within the first magnetic 
Brillouin zone given by | A;^ | < tv/Q and \k y \ < tt/2. 

Thus we obtain the eigenvalue equation, M^P = E^f, 
where ^ = (ipi t 1,^1, 2,^2, 1,^2, 2, • • • , ipQ, i,ipQ, i) T and 
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We calculate the energy eigenvalues for all the values of 
k in the first magnetic Brillouin zone, \k x \ < tt/Q and 
\k y \ < 7r/2, by directly diagonalizing the 2Q x 2Q Hamil- 
tonian matrix M(fc). As indicated in Fig. 1, the system 
has the checkerborad translational symmetry. That is, 
the system is invariant under the lattice translation by 
two lattice constants along either the x or the y direc- 
tions, or under the translation (n, m) — * (n. + 1, m + 1) 
along the diagonal. Thus one expects that, under the 
above transformations, the energy spectrum obtained 
by the eigenvalue problem with the Hamiltonian matrix 
M(fc) should remain the same, which is not obvious as 
seen from Eq. (0. In Appendix A, we prove that the 
Hamiltonian matrices before and after the translations 
are identical up to a shift in k y and thus give the same 
energy spectra. 

We show that there are several general symmetries of 
the spectra in the E — diagram, which can be used 
to reduce the amount of calculations. Similar discus- 
sion for the systems without field modulations can be 
found in Ref . |7J . In the following, the collecton of en- 
ergy subbands at a flux (per plaquette) modulated by 
(A^, Ay) is denoted by E((f>, A^,, Ay). It has the follow- 
ing symmetries^ 

I. E{f, A , Ay) = E(-cj), -A*, A y ) 

This follows from using two (three-dimensional) co- 
ordinate systems which are mirror images of each 
other with respect to the x — y plane. The physics, 
and hence the energy spectra, should be the same 
in these two frames with opposite handnesses. 

II.E(0,A ,Ay) £(-0,A ,Ay) 
£(-0,-A ,-Ay) 

The first equality follows from rotating the (three- 
dimensional) coordinate frame by 180 degrees 
around either the x-axis or the y-axis; the second 
is from shifting the coordinate by one plaquette 
along either the a;-axis or the y-axis. 

III. E(4>, A , Ay) = E{4> + 20 o , A , Ay) 

This results from the following two facts: (i) the 
smallest hopping loop for electrons encloses half 
of a plaquette; (ii) the Aharonov-Bohm phase for 
this closed loop is unchanged after adding one flux 
quantum to within this loop. 

IV. E(<j>, A^, Ay) = ^.E{4> + O , -A , Ay), where the 
operator 3? flips the spectrum with respect to the 



horizontal E — line. 

This follows from the two transformations: (i) 

fn,m -> (-l)" +m /n,m and -> + O ; (ii) f n ,m -> 

/n,m+i and A^ — > — A^. It can be shown that the 
overall sign of the Hamiltonian in Eq. changes 
after this two transformations, thus the spectra 
should have symmetry IV after shifting by one 
0o and reversing the direction of A^ . 

V. E(4>, A , Ay) = E{4> - 0o, A - 1, Ay) 

This is because of the freedom in shifting the fluxes, 
(0A, 0s) -> {(j>A -20o, 0s), where cj) A = (f> + 2-irA lj , 
and 0_b = — 2-7rA0 are the fluxes through each 
plaquette of the A and B sublattices respectively. 

By combining symmetries I and II, it is not difficult to 
see that E((p, A^,Ay) should remain unchanged when 
the sign of any of its arguments, 0, A^, or Ay, is 
changed. From symmetries III and II, we have -B(0o + 

0, A*, Ay) = E(<j> - 0o, A*, Ay) = £(0o - 0, A , Ay). 

That is, the distribution of E((f) in the E — cf> diagram 
has a mirror symmetry with respect to the vertical line 

= O . 

Finally, we show that, in the E — diagram, it is suf- 
ficient to plot the spectra within the range < < 0o/2 
only. The reason is as follows: from II, III, IV, and 
the freedom to flip the signs of A^, and Ay without 
changing the spectrum, it is clear that, for fixed val- 
ues of A^ and Ay, it suffices knowing the spectrum 
within the interval [0, 0o). Moreover, from IV and I, we 
have £(0 o /2 + 0, A , Ay) = 5i£(0 - o /2, -A , Ay) = 
5fti?(0o/2 — 0, A^, Ay). Therefore, the spectrum along 
the whole flux-coordinate can be determined by the E((f>) 
within the interval [O,0o/2). 

III. MAIN FEATURES OF THE SPECTRA 

In the discussion below, all energies are in units of 
t\. Besides ii(=l), there are three adjustable parameters 
in the present generalized Hofstadter model on a square 
lattice: A^, Ay, and t 2 - It is impossible to show all 
the results from the whole three-dimensional parameter 
spaced Therefore, we selectively report on certain sets 
of parameters with representative features. Notice that 
m Refs. PITliU - neither electric modulation nor NNN 
hoppings has been considered and the parameter space 
is one dimensional only. 
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First, we consider the effect of the checkerboard mod- 
ulation on the systems without NNN couplings t 2 , which 
have several symmetries in addition to the symmetries 
I~V listed above. 

EH. E(4>, A , Ay) = E(4> + O , A , Ay) 

When there is only NN hoppings, the period of the 
spectrum is one flux quantum since the smallest 
loop of hopping now is one plaquette, instead of 
half of the plaquette. 

EI. E(<f>, A , A v ) = ME{<f>, A , A v ) 

This results from symmetries IV and III', followed 
by flipping the sign of A^, which would not change 
the spectrum. Thus, E{<f>, A^,Ay) is symmetric 
with respect to the horizontal E = line when t 2 

EI. E(<j>, A , A y ) = E(cf> - cf> /2, A - 1/2, Ay) 

The argument is similar to the one leading to V, ex- 
cept that now <pA can be shifted by one 4>o without 
altering the Aharonov-Bohm phase of a closed-loop 
hopping. 

In Fig. 2, the spectrum for a checkerboard modulation 
with (A<£,Ay) = (0.1,0) is presented. The spectrum is 
indeed symmetric with respect to the E = line, ac- 
cording to symmetry IV'. Furthermore, because of the 
symmetry V and the freedom to flip the signs of the argu- 
ments, wehave£ , (0 o /4 + 0, A^, Ay) = E(<ft-<f>o/4,l/2- 
A , Ay) = £(0 o /4 - <j>,A<f> ~ 1/2, Ay). Therefore, af- 
ter being reflected by the vertical line at 0o/4, Fig. 2 
with A^ = 0.1 is identical to the Hofstadter spectrum 
for A = 0.4.21 

In Fig. 3, a checkerboard modulation with (A^, Ay) = 
(0.1,0.1) is considered. Without NNN hoppings, this 
Hofstadter spectrum retains the same symmetries as in 
Fig. 2. However, a distinctive ^-independent energy gap 
with a magnitude E g = 2Ay appears in the middle [also 
see Fig. 6(a)]. This is true with or without adding the 
modulation A^. First, it is not difficult to understand 
why the spectrum splits to two groups in energy: they 
originate from the two Bloch bands at cf> — due to the 
checkerboard modulation of the scalar potential. What 
is surprising is that the magnitude of the gap remains 
a constant for different </>'s and A^'s. It is no longer 
a constant as long as NNN couplings are included (see 
Fig. 4). A proof of the existence of the flux-independent 
gap is given in Appendix B, where it is shown that 
there exists a very simple relation between the spectra 
with and without electrostatic modulation Ay. That 
is, E(4>, A , Ay) = ±[E{(j), A^O) 2 + Ay] 1 / 2 . It can be 
checked that the spectra in Fig. 2 and Fig. 3 do obey 
this relation in details. 

The constancy of the energy gap in the limit of small 
flux 4> can be understood in the following semiclassical 
pictured The energy bands with vanishing widths as 
<j> — > in Fig. 3 are the cyclotron energy levels of the two 
parent bands at <fi — 0, which have the energy dispersions 



E±(k) = ±[(cosfca; + cosfc a ) 2 + A 2 ,] 1 / 2 if A^ = 0. It can 
be shown that, near the two inner band edges with ener- 
gies E± = ±Ay, the cyclotron effective masses approach 
infinity. Therefore, the position of the lowest Landau 
level approaches the lowest possible energy at the band 
edge and does not depend on the uniform magnetic field. 

When NNN couplings are included, the spectrum im- 
mediately lose the mirror symmetry with respect to the 
horizontal E — linei£ If only two of the three parame- 
ters are nonzero, then the spectrum remains fractal but 
distorted. When all three parameters, A^Ay, and t%, 
are nonzero, the subbands become significantly wider in 
most, but not all, of the regions. A typical example is 
shown in Fig. 4. The extent of widening varies as the 
parameters are varied. Because of the widening, the elec- 
trons are more delocalized, and become more mobile in 
transport. 

There is a surprising exception to the asymmetry re- 
sulted from NNN hoppings: the symmetry is restored 
again when A^ = 0.5, even if both Av and t 2 are nonzero. 
For example, the symmetric spectrum shown in Fig. 5 is 
for (A 4> ,Av,t 2 ) = (0.5,0.3,0.7). The existence of such 
a symmetry can be proved as follows. From the sym- 
metries III, IV, and V listed above, and the freedom to 
flip the signs of the arguments, it can be shown that 
E(<f>, A^, Ay) = $tE(4>, 1 — A$, Ay), which is a far less 
apparent symmetry since it relates two systems with dif- 
ferent strengths of flux modulation. It is clear that when 
A^ = 0.5, the spectrum has to be symmetric with respect 
to the line E = 0. 

In Fig. 6, we demonstrate how the continuous varia- 
tion of Ay and A^, influence the spectrum. The uniform 
flux and the NNN couplings are fixed at the values of 
p/q = 2/5 and t 2 = 0. In principle, there should be 
Q = 2q = 10 bands at this value of the flux. However, 
in Fig. 6(a) with (Ay, A^) = (x,0), where x G [0,1] 
is the value of the ^-coordinate, only 6 bands are ob- 
served. In fact, each of the upper two and lower two 
bands is itself formed by two overlapping subbands. We 
can also see that the band gap in the middle is indeed 
proportional to Ay, as mentioned earlier. In Fig. 6(b), 
(Ay, A^) = (0,x), where x e [0,1] is again the value of 
the x-coordinate. There is almost no similarity between 
(a) and (b). The overlapped subbands in Fig. 6(a) are 
split by a nonzero A^ and become very thin in most of 
the regions. On the other hand, the band in the mid- 
dle is thick and is actually composed of two subbands. 
The increase of flux modulation also induces many band 
crossings. In addition, there is an apparent symmetry 
£(0,A^,Ay) = E(((>, 1 - A , Ay). In Fig. 6(c), both 
Ay and A^ are nonzero and have the same numerical 
value. It has mixing features from (a) and (b), but the 
magnitude of the energy gap in the middle is not altered 
[comparing with (a)] by the nonzero A^. Such a contin- 
uous tuning of the band structure might be realized in 
the future using the optical lattices formed by quantum 
optical meansi^i 
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IV. SUMMARY 

The studies of Hofstadtcr spectrum have evolved from 
pure academic curiosities to accessible experimental in- 
vestigations. It is a basic physics problem involving sim- 
ple interplay between a lattice and a magnetic field. Be- 
cause of its general setting, it is not surprising to find 
counterpart problems in different physical systems, such 
as the quantum Hall system, the type-II superconductiv- 
ity, and the two-dimensional antiferromagnetism. Moti- 
vated by a study on the frustrated antiferromagnetism, 
and the recent experimental advances, we study the Hof- 
stadter problem with checkerboard modulations in de- 
tails. In this paper, the spectra are found to have sev- 
eral flux-related symmetries with respect to the change 
of <j) and A<£. One unanticipated symmetry occurs when 
A$ = 1/2. At that value, the spectrum are symmet- 
ric with respect to the E = line even in the presence of 
NNN hoppings. In the absence of NNN hoppings, we find 
a flux-independent energy gap induced by electric mod- 
ulations. Furthermore, a simple connection between the 
spectra for bipartite systems with and without electric 
modulation is discovered. More detailed aspects of the 
spectra are not investigated in this paper, however, such 
as the change of the fractal measures in the — Ay — £2 
parameter space. Such a study would reveal different 
phases in this space, as was done by Han and coworkers 
on the systems in a uniform magnetic fields The most 
general problem, when the superlattices of modulation 
can have the symmetries of their own, is considerally 
more involved. This study offers a starting point for re- 
searches in this direction. 
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APPENDIX A: PROOF OF THE 
CHECKERBOARD-TRANSLATION SYMMETRY 
OF THE SPECTRUM 



cal to the original one with another value of k y , i.e., k y — > 
k y + 2(f). Thus the whole energy spectrum within the 
first magnetic Brillouin zone remains the same. Third, 
the matrix elements in the Hamiltonian matrix can be 
shown to obey the following identities: a x A n+ i(k y )a x — 
A„(fc y + 4>) and a x B n+1 (k y )cr x = B n (k y + 0), where a x 
is the Pauli matrix. By using these identities, one can 
prove that, under the lattice translation (n, m) — > (77 + 
1, m + 1) such that ip n ,m(k x ,k y ) — > ip n +i, m+i{kx, ky), 
the new Hamiltonian matrix again becomes identical to 
the original one with a shift k y — > k y + <f>. Hence we con- 
clude that the energy spectrum is indeed invariant under 
the checkerborad translation. 



APPENDIX B: PROOF OF THE EXISTENCE OF 
THE FLUX-INDEPENDENT ENERGY GAP 

For ti = 0, our model is a nearest-neighbor-hopping 
model on a bipartite lattice. Therefore, we can rewrite 
the Hamiltonian in Eq. (Q as 



# = ({/!}, {/!» 



A V I 
2?t 



V 

-Ay/ 



Ua} 



(Bl) 



where I denotes the identity matrix, {/a} = {fn.m \ n + 
to is even} is a set of fcrmion operators for sublattice A 
and {/b} = {f n ,m I n + m is odd} is for sublattice B. 
When Ay = 0, the Schrodinger equation is 



V\(*a 
2>t U B 



E 



(B2) 



where Eg is the eigenvalue for the system with Ay = 0, 
and ($>a,&b) T is the corresponding eigenvector. From 
them we can construct the eigenstates for the original 
problem: 



(B3) 
(B4) 



with the corresponding eigenvalues E± = ±y/E^ + A 
because 



2 

v 



In this appendix, we show that the energy spectrum 
obtained from Eq. JUJl does respect the checkerboard 
translation symmetry. 

First, because there is no m-dependence of the ma- 
trix elements in Eq. @, the Hamiltonian matrix and 
therefore the spectrum are unchanged under the lat- 
tice translation 777 — > 777 -|- 2 such that ipn, m{k X7 ky) — ► 
■0 n , m +2{k x , k y ) . Second, from Eqs. (J5J and J§J, one can 
show that the matrix elements in the Hamiltonian ma- 
trix satisfy the relations A n+ 2(fcj,) = A n (k y + 2<p) and 
B„+2(fcy) = B„(fc y + 2(f)). Therefore, under the lat- 
tice translation n — > 77 + 2 such that ip n< m (k x , k y ) — » 
ipn+2, m(k x , k y ), the new Hamiltonian matrix for the 
eigenvalue problem after transformation becomes identi- 



AyJ V 



*± =±JE^ + A^ ± . 



(B5) 



Therefore, the energy spectrum is symmetric with respect 
to the horizontal E = line as mentioned in Sec. III. The 
positive-energy and the negative-energy parts are sepa- 
rated by an energy gap 2Tj\E \l lin + A y , where |S |min 
is the minimum value of |22o| at given <f> and A^. Since 
it has been shown that zero-energy modes exist in the 
absence of Ay for all flux values we have |i?o|min = 
for all values of (f> and A^. Consequently, the magnitude 
of the energy gap in the presence of Ay should be 2 Ay, 
independent of the values of <fi and A^. 
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FIG. 1: A square lattice with checkboard field modula- FIG. 3: The Hofstadter spectrum with the following param- 

tion. The magnetic flux through the left (right) plaquette eters: (A^, Ay, £2) = (0.1,0.1,0). 

is <f> + 2-kA,/, (4> — 2-kA^,). The scalar potentials at the lattice 

points indicated by solid and empty dots are Ay and —Av 

respectively. 



FIG. 2: The Hofstadter spectrum with the following param- 
eters: (A^Ay.fa) = (0.1,0,0). 



FIG. 4: The Hofstadter spectrum with the following param- 
eters: (A , Ay, t 2 ) = (0.1,0.1,0.5). 
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FIG. 5: The Hofstadter spectrum with the following param- 
eters: (A^,A v ,t 2 ) = (0.5,0.3, 0.7). 



FIG. 6: The variation of band widths as (a) Ay, (b) A^>, and 
(c) both Ay and A^> are tuned. The values of 4>/2n — 2/5 
and ti — are fixed. 
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